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We propose a new method for atomic-scale imaging of spatial current patterns in nanoscopic 
quantum networks by using scanning tunneling microscopy (STM). By measuring the current flowing 
from the STM tip into one of the leads attached to the network as a function of tip position, one 
obtains an atomically resolved spatial image of "current riverbeds" whose spatial structure reflects 
the coherent flow of electrons out of equilibrium. We show that this method can be successfully 
applied in variety of network topologies, and is robust against dephasing effects. 
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How does a current flow through a nanoscopic system? 
The answer to this question is not only of fundamental 
interest for our understanding of quantum phenomena 
at the nanoscale, but also of great importance for the 
further development of nano-electronics and the contin- 
ued miniaturization of electronic devices. In mesoscopic 
systems, such as quantum point contacts [IH3], quan- 
tum rings [4] and DNA [5], this question has been suc- 
cessfully answered by imaging spatial current paths us- 
ing a scanning probe microscope (SPM) [THZj. However, 
in nanoscopic systems with sizes of tens of nanometers, 
this question has only been addressed theoretically so 
far [8 -10 . The main experimental challenge here arises 
from the SPM's perturbing electrostatic potential, which 
yields a spatial resolution that is insufficient to detect the 
atomic-scale varying current patterns predicted to exist 
in nanoscopic systems [8UT0]. Developing an experimen- 
tal method that allows atomically resolved imaging of 
spatial current paths is therefore of great importance for 
understanding charge transport at the nanoscale, open- 
ing new possibilities to explore quantum Hall physics [11] 
or weak localization effects [12] in nanoscopic systems. 

In this Letter, we propose such a method by demon- 
strating that atomic-scale imaging of spatial current 
paths in nanoscopic quantum networks can be achieved 
using scanning tunneling microscopy (STM) [13]. In par- 
ticular, by measuring the current flowing from an STM 
tip into one of the leads attached to the network as a 
function of tip position, one reveals atomically resolved 
"current riverbeds" whose spatial structure reflects the 
coherent flow of electrons out of equilibrium. We show 
that this method can be successfully applied in a vari- 
ety of network topologies, and reproduces even complex 
current patterns arising from the presence of constric- 
tions. It is also robust against dephasing effects, pro- 
viding correct current images even in systems where the 
mean- free path is only a few lattice constants. Finally, we 
demonstrate that the form of the network's conductance 
is an experimentally accessible criterion for the success of 
the imaging method. Thus, the proposed STM imaging 
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FIG. 1. (color online) Position of the chemical potentials of 
the leads, the STM tip, and the network for cases I and II. 



method provides a promising new approach for exploring 
charge transport in nanoscopic systems. 

We study the spatial current patterns in a nanoscale 
quantum network that is attached to two leads and de- 
scribed by the Hamiltonian 
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Here, c\ a and d\ a create an electron with spin a at site 
r in the network and leads, respectively, with t, t\ and th 
being the respective hopping elements. While the net- 
work sites can in general represent atoms, molecules or 
more complex structures, we assume here that each site 
contains only a single relevant electronic level. The last 
two terms describe a local electron-phonon interaction in 
the network with a phonon mode of energy cuq. 

For different chemical potentials, [1l,r in the left and 
right leads (see Fig. [I] case I), a spatial current pattern 
is obtained by computing the current, I r y, between ad- 
jacent sites i*,!* 7 in the network via [T4| [T5] 
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where, G K is the full Keldysh Green's function matrix of 
the entire system. For a non-interacting (g = 0) network, 

one has G K = (l - gH)' 1 g K (l - tg -)' 1 with g K (u) = 
2i [1 — 2uf(w)] Im [g r (uo)\. Here, t is the symmetric hop- 
ping matrix and the matrices g r ^ a ' K contain the decou- 
pled (t, th = 0) retarded, advanced and Keldysh Green's 
functions of the network sites [with g r = 1/(lj + i5), S = 
+ ] and leads. The latter are computed using a renormal- 
ization procedure [16 . hp contains the Fermi functions 
of the network and leads. Below, we use for concrete- 
ness t = t h = t u S = 0.0002t, A/i = fi L ~ Vr = 10~H, 
and k B T = 10~ 5 t. In the limit AV = Afi/e and 
temperature T = 0, one obtains from Eq.Q [TO] 
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where /i c = (ml + Vr)/2 and iVx, is the density of states 
(DOS) in the left lead. 

To image the spatial current patterns obtained from 
Eq.Q, we compute the current, Ir^r(T), flowing from 
an STM tip (held at potential jut) into the left (L) or 
right (R) lead as a function of tip position T (see Fig. [I] 
case II). Here, Ht = —tr E r a { c t,crfo- + H.c.) describes 
the electron tunneling between the STM tip and a single 
network site with ft creating an electron in the tip, and 
one appropriately extends G K in Eq.fji]). To probe the 
same states participating in the charge transport as in 
case I, we set the bias between the leads to zero, i.e., 
Hr,r = Mo, and require AVr = (/jlt — Mo)/e = AV and 
ji^ = (/io + Mt)/2 = \i c with AV and \i c from case I 
(see Fig. [TJ. Assuming, for example, that the left lead is 
connected to a single network site L only, one obtains for 
AV T -> and T = 



I L (T) = 2—t 2 T t 2 h N L N T \G r T:Ij \ 
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where N T is the DOS in the STM tip. Note that in 
contrast to i^, the total tip current, I t i p = II + -Tr 5 fo r 
AVt —> yields a spatial image of the network's density 
of states at uj = jj^ . In the experimentally realized weak 
tunneling limit [13] , the STM only probes but essentially 
does not perturb the network's electronic structure, and 
is thus less intrusive than the SPM imaging method [T]-[7] . 

We begin by considering the spatial current pattern, 
I r y, obtained from Eq.Q (case I) in a non-interacting, 
semi-infinite network with N y = 15 rows that is con- 
nected to a narrow left lead through a single site, L (see 
Fig. [2]). For \i c = [Fig. [2|a)], the Fermi velocity of 
the states involved in the current transport points along 
the lattice diagonal, resulting in a diagonal current path 
that bounces off the sides of the quantum network [TU] . 
A comparison of this current pattern with II obtained 
from Eq.Q (case II) shown in Fig. |2|b) demonstrates 
that II provides an atomically resolved image of the cur- 
rent paths in the network, a result which also holds when 
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FIG. 2. (color online) (a) Spatial current pattern, I r y, and 
(b) I L (T) for fi c = 0. (c) I r y, and (d) I L (T) for /x c '= 2 At. 
All currents are normalized. 



the network is gated, and fi c is shifted, as exemplified in 
Figs. [2|c) and (d) where we present the spatial current 
pattern and II, respectively, for \i c = 2 At. 

A physical understanding of the good agreement be- 
tween I Y y and II can be gained by considering how the 
current flows from the STM tip through the network into 
the left lead. Consider, for example, the case when the 
STM tip is located above a network site that is part of 
the coherent "current riverbed" of Fig. [2^ a). In this case, 
as shown in Fig.[3ja), the current injected from the STM 
tip utilizes the current riverbed, resulting in a large Ir 
flowing from the tip to the left lead. In contrast, when 
the tip is positioned away from the riverbed, as shown 
in Fig.|3jb), the current injected by the STM tip cannot 
flow to the left lead, but only to the wide right lead (at 
infinity), since only the latter can be connected to the tip 
location by a diagonal path. While the spatial form of 
current riverbeds varies with the position of the leads and 
the gating of the network, this physical explanation for 
the success of the imaging method holds for all cases we 
have considered. A more technical understanding of the 
good agreement between I r y and Il : is obtained by not- 
ing that their analytical expressions in Eqs.([3| and Q 
are determined by (different) combinations of the real 
and imaginary parts of the non-local Green's function, 
G^l- A spatial plot of these two quantities shown in 
Figs. [3fc) and (d) [for the case of Fig. [2|a)] reveals that 
their spatial structures are similar on the atomic scale. 
As a result, Ir [see Eq.Q] is strongly peaked in those 
regions where I r y ^ (the latter being predominantly 




FIG. 3. (color online) (a),(b) Spatial pattern of current flow- 
ing from the STM tip to the leads for two different tip po- 
sitions (indicated by yellow triangles) and \±t — 5 x 10 _5 t, 
He = 0. (c) ReCr,L an d (d)ImG> jL f° r the case °f Fig.^a). 



determined by ReG^ L [10]), thus explaining the good 
spatial agreement between II and I r y- The same con- 
clusion also holds for the case of [i c = 2 At in Figs. [2|c) 
and (d). As follows from Eqs.Q and Q, and confirmed 
by an extensive survey of networks, a spatial structure of 
ReG^ L and ImG^ L that is similar on the atomic scale 
is a sufficient condition for II to provide an atomically 
resolved image of current paths (see discussion below). 

The STM method also images more complex current 
patterns, which occur, for example, in systems with mul- 
tiple constrictions, as shown for an infinitely long wire 
with N y = 13 rows in Fig. ga). While in the left and 
right regions of this network, the current follows a di- 
agonal path, a more complex current pattern exhibiting 
circulating current loops [10] exists in the central region. 
As shown in Fig. Qb), 1 L correctly images the current 
pattern in the central and right regions, but reflects the 
spatial form of the density of states at ji^ in the left re- 
gion. Conversely, Ir (not shown) correctly images the 
spatial current pattern only in the left and central re- 
gions. To understand why Il^r reproduce the current 
pattern only in parts of the system (while together, I l ,r 
yield a complete image of the current paths), we note 
that when the tip is located to the right of constriction 
1 [Fig. |4jb)], the current flowing from the tip into the 
left lead needs to pass through the narrow opening in 
constriction 1, whose role is therefore similar to that of 
the narrow left lead in Fig. [5|a). As a result, I L is large 
whenever the tip is located above the current riverbed 
of Fig. |4ja). In contrast, in the region left of constric- 
tion 1, the tip, independent of its location, can always be 
connected to the wide left lead (at infinity) by a diago- 
nal path, and I L ~ hip therefore images the density of 
states. We thus conclude that in order for II to correctly 
image the current pattern in a given region, the number 
of allowed current channels need to be restricted, either 
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FIG. 4. (color online) (a) Spatial current pattern and (b) 
I L (T) for fjic = 0. (c) Jl(T) for t L = 0.2. The scales in (b) 
and (c) are the same. Constrictions are shown in gray, (d) 
G(/jl c ) for several tL- 



through the presence of a constriction, as in Fig.|4jb), or 
of a narrow lead, as in Fig. [2|c). Similar conclusions also 
hold for more complex network topologies, such as the 
infinitely long cylinder of circumference N y = 54 shown 
in Fig. [5ja) possessing two constrictions rotated by a rel- 
ative angle of tt around the cylinder axis. Here, II [see 
Fig. [5^b)] correctly images the winding current pattern 
only in the central and right cylinder regions, while Ir 
provides the complementary current image. 

To make the imaging method a versatile tool in the 
study of nanoscopic networks, it is necessary to identify 
an experimentally verifiable criterion for the agreement 
(or the lack thereof) between II and the current pat- 
tern. Such a criterion is provided by the network's over- 
all conductance, G(/i c ), as can be exemplified by varying 
the hopping, t^, between sites in constrictions 1 and 2 
[see yellow dots in Fig. [4^c)] and the central region of 
the network shown in Fig. [4j As discussed above, for 




FIG. 5. (color online) (a) Spatial current pattern for /x c = 0, 
and (b) Il(T) for a cylinder (cut along the axis and flattened 
out). Constrictions are shown in gray. 



tL = t, the STM method provides a correct spatial im- 
age of the current pattern while at the same time, the 
conductance around \i c = [see arrow in Fig. (4^d)] is 
close to the (maximal) quantum of conductance, vary- 
ing only weakly with \i c . However, with decreasing t^, 
the agreement between the STM image and the actual 
current pattern (the latter remaining qualitatively un- 
changed) worsens in the central region and is completely 
lost around th — 0.2t [see Fig. [4^c)]. This break down 
is accompanied by a qualitative change in G(/i c ): it be- 
comes sharply peaked around fi c = 0, and its overall mag- 
nitude decreases. Both changes possess the same origin: 
for II —± 0, only a single state contributes to the charge 
transport at /i c = in the central region, implying that 
G(/i c ) becomes sharply peaked while the overall magni- 
tude of ImG^ L (/i c ) scaling with £ _1 becomes much larger 
than that of ReG£ l (aO- As a result, II is dominated by 
ImG and therefore does not correctly image the current 
pattern, whose spatial form is in general determined by 
ReG [TO] . We find that in general, the onset for the break 
down of the imaging method occurs when the widths of 
the peaks in G(/i c ) become less than the separation be- 
tween peaks [see = 0.5t in Fig. J4^d)]. In contrast, 
even for = 0.2t, II still correctly images the current 
pattern in the right region since due to its semi-infinite 
nature, ImG does not scale with £ _1 for — ^ 0, and 
the overall magnitude of ImG and ReG remain the same 
(though the overall magnitude of II rapidly decreases 
for tL < 0.2t). This relation between the success of the 
imaging method and G(/i c ) holds for all networks we have 
considered: whenever G(/i c ) is close to the maximal al- 
lowed conductance, and varies weakly around a given /i c , 
the STM method yields a correct current image in the 
entire network. On the other hand, if G(ja c ) is sharply 
peaked (or reduced from its maximal value due to dephas- 
ing, as discussed below), the agreement breaks down in 
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FIG. 6. (color online) Spatial current pattern and II for 
(a),(b) 7 = 0.002£ 2 and (c),(d) 7 = 0.032t 2 with n L = -fi R = 
0.03£, k B T = 5 x 10" 6 t = 10cj , and S = 10" 16 t. 



at least part of the network. Thus, G(/i c ) provides an 
experimentally verifiable criterion for the success of the 
imaging technique. 

To investigate how the ability of Il,r to image current 
patterns is affected by interactions and the resulting de- 
phasing, we consider the effects of the electron-phonon 
interaction in Eq.0[17j[18]. The full solution of the self- 
consistent Dyson-equation for G K,r for general g, ujq and 
temperature T is currently beyond our computational 
abilities, and thus reserved for future work. However, 
in the high-temperature approximation k B T ^> uo (with 
uo — » 0) [19 , the fermionic self-energy greatly simpli- 
fies to S K ' r = 7 G K ' r , where 7 = g 2 [1 + 2n B (wo)], and 
the Dyson equation can be numerically solved even for 
larger networks [10| . In Fig. [6] we plot the spatial cur- 
rent pattern and 4 for a (17 x 5) network attached to 
a narrow and a wide lead, and two different values of 
the dephasing parameter 7. In general, with increasing 
7, the current pattern evolves from that of the ballis- 
tic limit to that of a classical resistor network [10 . For 
small 7 up to 7 ~ 0.002t 2 [corresponding to a mean-free 
path of I « 13ao, where ao is the lattice constant], with 
conductance G w 0.99 (2e 2 /h), I L [Fig. ^b)] correctly 
images I Y Y > [Fig. |6|a)] which still largely resembles the 
current pattern of the quantum (7 = 0) limit. For larger 
7 = 0.032£ 2 with I « 3.3a and G « 0.93 (2e 2 /h), both 
ir )r / [Fig. [6|c)] and II [Fig. [6jd)] become spatially more 
diffuse, however, iz,(T) still reproduces the main features 
of the current pattern even for such a small mean-free 
path. This demonstrates that the STM imaging tech- 
nique is rather robust against dephasing effects even for 
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mean- free paths as small as a few lattice spacings. In- 
creasing 7 even further results in an overall suppression 
of G(/x c ), and an increasing disagreement between II and 
7 r r /, in agreement with the criterion formulated above. 
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